Some Jump Processes in Quantum Field Theory 

Roderich Tumulka* and Hans-Otto Georgia 



Abstract 

A jump process for the positions of interacting quantum particles on a lattice, 
with time-dependent transition rates governed by the state vector, was first con- 
sidered by J.S. Bell. We review this process and its continuum variants involving 
"minimal" jump rates, describing particles as they get created, move, and get 
annihilated. In particular, we sketch a recent proof of global existence of Bell's 
process. As an outlook, we suggest how methods of this proof could be applied 
to similar global existence questions, and underline the particular usefulness of 
minimal jump rates on manifolds with boundaries. 

MSC (2000). 60J75 , 81T25 

1 Introduction 

This contribution deals with Markov jump processes Qt describing the positional time 
evolution of finitely many interacting quantum particles. These processes are char- 
acterized by a specific form of time-dependent jump rates induced by the Schrodinger 
equation for the quantum state vector \l/ t of the underlying quantum field theory (QFT). 
As a typical example, suppose that the particles live in the physical three-space M 3 . The 
process Qt then takes values in the space of all finite subsets of M 3 , the configuration 
space of a variable number of identical particles (corresponding to the Fock space of 
QFT), and the jumps of Q t describe the creation or annihilation of particles; between 
these stochastic jumps, Qt evolves deterministically according to some ordinary differ- 
ential equation governed by tyf Alternatively, one may think of quantum particles on 
a lattice; the jumps of Qt then record all changes of the particle configuration. We will 
portray several processes of this type, present some common principles, and in particu- 
lar discuss some results of two recent papers |GT03t rDGTZ03c] . the work on which was 
supported by the DFG Priority Program. 

As the state vector \l/ t determining the jump rates follows the time-dependent Schro- 
dinger equation, the jump rates themselves are explicitly time-dependent, so that the 
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processes Q t considered here do not admit an invariant measure. However, the jump 
rates are designed in such a way that Q t does admit an equivariant measure, namely the 
quantum distribution l^l 2 , which means that Q t has distribution \^t\ 2 at & n y time t. 
This is the key feature justifying the particular form of the jump rates, and on the other 
hand the main fact on which one can build the analysis of these processes. So, the issue 
here is not the analysis of distributional properties of a given process, but the converse: 
the equivariant distribution is given, and the objective is to prove the existence of the 
associated process, and to check that it really does have the equivariant distribution. 
In |GT03j . we carried out this program for the case of a discrete configuration space, 
including in particular a lattice model proposed by J.S. Bell Bel86]; the main arguments 
will be sketched in Sect. 4. 

From the probabilistic viewpoint, one has to overcome two difficulties. First, the 
transition rates exhibit singularities, in that they become ill-defined at certain time 
and space points. One has to show that the process avoids these singularities. The 
second (and more important) task is to rule out the possibility of explosion, i.e., the 
accumulation of infinitely many jumps in finite time. Due to the unbounded growth 
of the rates near the singularities, the standard methods fail, and one has to use the 
particular relation between transition rates and equivariant distribution. 

Besides the results on the discrete case mentioned above, we will also describe the 
continuum analogues of Bell's process investigated in |DGTZ03a] IDGTZ03cj : as a spe- 
cial case these include Bohmian mechanics Boh52t IBel66| IBDDG Z951 IDurOlj for the 
appropriate Hamiltonian with a conserved number of particles. On the basis of what 
we learned from our existence proof for Bell's model, we also propose here some new 
methods for proving the existence of Bohmian mechanics. 

Let us now discuss how the models considered here relate to the topics of other 
articles in this volume. First, the existence problem for a model of quantum field theory 
is also the subject of the contribution of S. Albeverio, Y. Kondratiev, M. Rockner, 
and T. Pasurek. The issue there is the existence, and uniqueness, of Euclidean Gibbs 
measures for infinitely many interacting quantum spins. These concern an equilibrium 
setting, and time appears only via path integration to make the connection with the 
quantum states. The difficulty there is the infinite number of spins, requiring particular 
assumptions on the interaction. By way of contrast, the models considered here involve 
only finitely many particles, but in a nonequilibrium situation, and we do not need any 
particular assumptions on the interaction. 

From the methodological side, there is a closer connection to the contributions deal- 
ing with population biology, in particular that by R. Hopfner and E. Locherbach. The 
similarities concern the creation and annihilation of particles in the vicinity of other 
particles, and the necessity of proving non-explosion. Also, the space-dependence of 
the reproduction rates in Hopfner and Locherbach's article implies non-exponential life- 
times of individuals, just as the time- dependent jump rates considered here imply non- 
exponential interjump times. However, in our case the paths between the jumps are 
smooth and deterministic. 

This note is organized as follows. In Sect. 2 we derive the fundamental formula 
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(JHJ) for minimal jump rates, defining the jump process associated with a certain type 
of Hamiltonian. This involves consideration of the equivariant probability distribution 
(J3J) and probability current (JHJ) provided by quantum theory. In Sect. 3 we explain 
the connection with Bohmian mechanics and with Bell's model. We also describe the 
process for a concrete example QFT, introduced in [E)GTZ03aj. In Sect. 4 we sketch the 
global existence proof for the discrete case, including Bell's model, that we developed in 
|GT03j . In Sect. 5 we point out how the methods of [GT03 could be adapted to other 
global existence problems. In Sect. 6 we indicate some perspectives for future research 
concerning a process for quantum theory on a manifold with boundaries, and the special 
role the minimal jump rate (JHJ) plays for this process. 



2 Jump Rates Induced by Schrodinger Equations 

We now introduce the class of jump processes we are concerned with, starting with a 
general framework. For our purposes, a quantum theory is abstractly given by a Hilbert 
space containing the state vectors, a one-parameter group U t of unitary operators 
on Jj? defining the time evolution 

% = Wo (1) 

of the state vector, and a measurable space of configurations describing the 

locations of particles. Q is tied to Jtf via a projection-valued measure (PVM) P(dq) 
on Q acting on Jt?, i.e., a mapping from the a-algebra & to the family of projection 
operators on ffl that is, like a measure, countably additive (in the sense of the weak 
operator topology) and normalized, in that P(Q) = I, the identity operator on Jff. If 
Jt? = L 2 (Q, C k ) with respect to some measure on Q then Jf? is equipped with a natural 
PVM, namely P(B) being multiplication by the indicator function of the set B. In 
nonrelativistic quantum mechanics, another way of saying this is that P is the PVM 
corresponding to the joint spectral decomposition of all position operators. 
By Stone's theorem, the unitary operators Ut are of the form 

U t = e' iHt/h (2) 

with H a self-adjoint operator on Jff, called the Hamiltonian. Equations (JTJ and (j2*|) 
together correspond to the formal Schrodinger equation 

ih-^ = H%. (3) 

We will show how this Schrodinger equation, together with the PVM P, gives rise to a 
natural Markov process on Q. 

In this section we focus on the case in which this Markov process is a pure jump 
process. (Roughly speaking, this will require that the Hamiltonian is an integral operator 
in the "position representation" defined by P; differential operators will be considered 
in Sect. 3.) So we ask: Is there any distinguished Markovian jump process (Q t ) on Q 
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describing the evolution of the particle configuration, and what are its transition rates? 
To answer this question we note that the quantum theoretical probability distribution 
of the configuration at time t is given by 

ir t {-) = (%\P{-)%). (4) 

(We generally assume that H^oll = !•) It is therefore natural to stipulate that 7t t is 
equivariant for (Qt), meaning that Qt has distribution Tr t at every time t. Can one 
choose some (time-dependent) transition rates (a t ) for (Q t ) to satisfy this requirement 
of equivariance? Yes indeed, in view of (j3J) the time derivative of ir t is given by 

n{B) = \lm{%\P{B)m t ) = f J t (B,dq'), (5) 

where 

J t (B, B') = I Im (%\P(B)HP(B')%) (6) 

is the quantum theoretical current between two sets B, B' E & . On the other hand, 
suppose (Q t ) is a pure jump process on Q jumping at time t with rate a t (B,q') from 
q' E Q to some configuration in B E JP. Then its distribution p t = P o Q^ 1 evolves 
according to the equation 

p t (B) = [ a t (B,q)p t (dq)- [ a t (Q, q) p t (dq). (7) 



To satisfy the condition of equivariance we need to find jump rates a t such that the 
right-hand sides of the evolution equations (0) and (JJJ) coincide whenever p t = n t . We 
see that this is the case when at is given by the Radon-Nikodym derivative 

Jt + (dq,dq') [llm(^t\P(dq)HP(dq')^t)} + 

a t (dq,q) - 



Tr t (dq>) (%\P(dq')%) 

of the positive part J t + of J t in its second variable q', provided this makes sense. For, 
the antisymmetry of Jt then implies that 

o- t (dq, q')nt(dq') - a t (dq', q)n t (dq) = J t (dq, dq') . 

To make formula (jHJ) meaningful one needs some assumptions which roughly require 
that H is an integral operator in the position representation given by P, and (Q, 
is standard Borel. This is discussed in detail in jt)GTZ03cJ, where (jEJ) was written 
down for the first time in this generality; special cases had been used before in [Bel86, 
Sud87, DGTZ03a^. For the precise conditions we refer to Theorem 1 (Sect. 4.1) and 
Corollaries 1-3 (Sect. 4.2) of |DGTZ03c| . Under these conditions, formula (jHJ) can (and 
has to) be read as follows: A priori, J t is a signed bi-measure on J^x (a measure in each 
of the two variables q, q'). This has to (and then can) be extended to a signed measure 
on the product a-algebra & ® & . The positive part in (jHJ) is then to be understood 
in the sense of the Hahn- Jordan decomposition of this extended measure. Next one 
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notes that, for each B G ^ , J t (B, ■) <C ir t because P(B')^/ t = whenever n t (B') = 0. 
One can thus form the Radon-Nikodym derivatives a t (B, ■) = Jf~(B, dq 1 ) / TT t (dq') , which 
finally have to be chosen in such a way that at becomes a measure kernel. 

According to our derivation above, the transition rates (JHJ) have been chosen to satisfy 
the requirement of equivariance. There was, however, still some freedom of choice. The 
particular rate (JHJ) is singled out by the following additional facts. 

1. Suppose there exists a jump process (Q t ) on Q with rates (JHJ). As is evident 
from the arguments above, the net probability current of (Q t ) between two sets 

B,B'e&, 

jt(B, B') = lim - (p(Q t G B', Q t+£ e B) - P(Q t G B, Q t+£ G B')) , 

then coincides with the quantum theoretical current Jt defined by (JHJ). Conversely, 
if (Qt) is any pure jump process having initial distribution tcq at time 0, some jump 
rates a t and probability current j t = J t , it turns out that necessarily & t (dq,q') > 
a t (dq, q') |RS90L IBD991 ID(;TZ03cj . This follows from the minimality of the Hahn- 
Jordan decomposition. The rates (JHJ) are therefore called the minimal jump rates, 
and a process with rates (JHJ) is distinguished among all processes with the right 
probability current by having the least frequent jumps, or the smallest amount of 
randomness. 

2. Always one of the transitions q' — > q or q — >• q' is forbidden. More precisely, for 
every time t there exists a set Sf G & ® & which, together with its transposition 
Sf, covers Q x Q (except possibly the diagonal), and such that 

a t({l '■ (q,q') £ Sf},q') = for 7r t -almost-every q. 

Indeed, by the anti-symmetry of J t , its positive and negative part J t + and J 4 ~ 
admit supports and that are transpositions of each other, whence the 
result follows. 

Put more simply, the mechanism is this: When the current J t (dq,dq') is positive, 
meaning that there should be a net flow from dq' to dq, then a t (dq,q') > and 
<7 t (dq',q) = 0, i.e., only jumps from q' to q are allowed; the converse holds in the 
case of a negative current. Under all rates with j t = Jt, the minimal rates (JHJ) are 
characterized by this property. 

3 Bohmian Mechanics and Bell- Type QFT 

In this section we discuss three particular instances in which jump rates of the form (JHJ) 
play a significant role. 

A. Bohmian mechanics as continuum limit of jump processes. Consider nonrelativis- 
tic quantum mechanics for N particles: the configuration space is Q = M. 3N , the Hilbert 
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space J? = L 2 (R 3N , C k ) and the Hamiltonian 



N h 2 

H = -Y,^.^ + V ^---^n) (9) 

i=i * 

with Aj the Laplacian acting on the variable Xj, the mass of the z-th particle, and 
V the potential function (possibly having values in the Hermitian k x k matrices). We 
obtain a Markov process on the configuration space in the following way: first discretize 
space, i.e., replace R 3 by a lattice A = eZ 3 and the Laplacian Aj by the corresponding 
lattice Laplacian Af. We then can consider the jump process Q\ on with rates (jSJ). 
As the lattice shrinks, e — > 0, one obtains |Sud87l IVin93j in the limit the deterministic 



process Q t satisfying the ordinary differential equation 

dQ M h VtVjVt 

ir = ^ M -^*r (g *' i '-'^ ) - (10) 

Here Q t) i is the i-th component of Q t , i.e., the position of the z-th particle, ^ t obeys the 
Schrodinger equation (jSJ) with Hamiltonian (jUJ), and $* $2 is the inner product in C fc . 
The process (jlUJl is known as Bohmian mechanics [Boh 52"l lBe!66| (BDDG Z951 IDurOlj . 
For a suitable other choice of jump rates |Vin93j . also satisfying j t = J t but greater than 
minimal, one obtains in the continuum limit e — > the diffusion process introduced by 
E. Nelson and known as stochastic mechanics |Nel85j . 

What makes Bohmian mechanics (or, for that matter, stochastic mechanics) partic- 
ularly interesting to quantum physicists is that in a Bohmian universe - one in which 
the particles move according to (jl(Jj) and the initial configuration is chosen according to 
the |\l/| 2 distribution - the inhabitants find all their observations in agreement with the 
probabilistic predictions of quantum mechanics - in sharp contrast with the traditional 
belief that it be impossible to explain the probabilities of quantum mechanics by any 
theory describing events objectively taking place in the outside world. 

B. Bell's jump process for lattice QFT. The study of jump processes with rates 
(|SJ) has been inspired by Bohmian mechanics, in particular by the wish for a theory 
similar to Bohmian mechanics covering quantum field theory. The first work in this 
direction was Bell's seminal paper [Bel86 . For simplicity, Bell replaces physical 3-space 
by a lattice A and considers a QFT on that lattice. A configuration is specified in his 
model by the number of fermions q(x) at every lattice site x. Thus, with the notation 
Z + = {0, 1, 2, . . .}, the configuration space is 



Q = r(A):={gGZ^:^g(x)<oo}, 



the space of all configurations of a variable (but finite) number of identical particles on 
the lattice. (In fact, he imposes a bound on the total number of particles and assumes 
that A is finite, but this is not necessary.) The PVM P( ■ ) that he suggests arises from 
the joint spectral decomposition of the fermion number operators N(x) for every lattice 
site, i.e., P(q) := P({q}) is the projection to the joint eigenspace of the (commuting) 
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operators N(x) for the eigenvalues q(x). The jump rate Bell uses is the appropriate 
special case of (JSJ): the rate of jumping from q' to q is 



For stu dies of Bell's process we refer to jSud87llVin931IBD99iran!?all(;ol03btll)(;TZ03cl 
GTQ3j, and for some numerical simulations and applications to [DR031 IDen03| . We 
return to it in more detail in the next section. 

C. Bohmian mechanics with variable number of particles. A third example of a 
process for a QFT was considered in (DGTZ03a1 . It arose from an attempt to include 



particle creation and annihilation into Bohmian mechanics by simply introducing the 
possibility that world lines of particles can begin and end. That is, the aim is to provide 
a generalization of the Bohmian motion (|lUj) to a configuration space of a variable 
number of particles. Here we describe this model in a simplified version. For the 
numerous similarities between our model process and Bell's discrete process, we called 
it a "Bell-type QFT." In DGTZ03q, methods are developed for obtaining a canonical 
Bell-type process for more or less any regularized QFT. 

A configuration of finitely many identical particles can be described by a finite count- 
ing measure on M. 3 . Since the coincidence configurations, those in which there are two or 
more particles at the same location, form a subset of codimension 3, they are basically 
irrelevant, and it will be convenient to exclude them from the configuration space. What 
remains, as the space of "simple configurations", is the set of all finite subsets of M 3 , 

r_,(M 3 ) = {gcK 3 :#(/< oo}. 

Under the physical conditions prevailing in everyday life, the most frequent type of 
particle creation and annihilation is the emission and absorption of photons by electrons. 
This can be described in a model QFT as follows. Particles (photons) move in a Bohmian 
way and can be emitted and absorbed by another kind of particles (electrons). For 
simplicity, we will assume here that the electrons remain at fixed locations, given by 
a finite set )j C I 3 ; the case of moving electrons is described in [UGTZ03a . The 
configuration space is thus the space of photon configurations, Q = r^(lR 3 ), and ^ t a 
square-integrable complex-valued function on Q; the Hilbert space J4? of these functions 
is known as the bosonic Fock space arising from L 2 (M 3 ). 

The Markov process Q t in Q has piecewise smooth paths. It obeys the deterministic 
motion (|1U|). interrupted by stochastic jumps. The process is piecewise deterministic in 
the sense that, conditional on the times of two subsequent jumps and the destination 
of the first, the path in between these jumps is deterministic. The jumps correspond 
to creation or annihilation of a photon near some point of 77; in particular, every jump 
changes the number of photons by one. The process is thus a special kind of a spatial 
birth-and-death process with moving individuals |Pre76j . 

The deterministic motion, during which the number of photons is kept constant, is 
defined by (jl(Jj) ; for simplicity, we deviate a little from the physical facts and assume that 
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the "photons" have a positive mass m p h- The rate for the transition q —>■ qUx := qU{x}, 
i.e., the creation of a new photon at the location x G M 3 \ q, has density (with respect 
to Lebesgue measure dx) 



I Im V*{q U x) (#<? + l) 1 ' 2 Even ^ x ~ V) **(?) 



a t (gUx,g) 



where <p : IR 3 — > R is a fixed function, a spherically symmetric, square-integrable po- 
tential supported by the ball of radius 5 > 0. Likewise, for any x G g the rate for the 
transition g — > q \ x := q \ {x}, i.e., annihilation of the photon at x, is 



These rates, together with vanishing rate for any other transition, are in fact a special 
case of (JHJ), for a suitable integral operator Hj in place of H. For the definition of Hj 
and the derivation of the above expressions from (|SJ) we refer to Sect. 3.12 of [DGTZ03cJ. 

Now, Hj is not the Hamiltonian of the relevant QFT, but its interaction part; i.e., 
the complete Hamiltonian is H = H + Hj, where H , the free Hamiltonian, is given by 



Observe that there is a correspondence between the splitting H = Ho + Hj and the two 
constituents of the process, the motion and the jump rates just given. Deterministic 
motion corresponds to Hq while the jumps correspond to Hi. Indeed, the minimal 
process, the one arising as a limiting case from jump processes with minimal rates 
(JHJ), associated with H alone is the continuous motion (jlUj) while the minimal process 
associated with Hj is the pure jump process with the above rates. 

This is an instance of the general rule of process additivity: If the minimal process 
associated with Hi has generator S£\^ t and the one associated with Hi has generator 
J2?2,* t , then the minimal process associated with H\ + H 2 has generator =Sfi,* t + -^2,* t > 
provided that the (formal) integral kernels of Hi and H 2 have disjoint supports. The 
sum of the generators of a deterministic flow and of a pure jump process generates the 
piecewise deterministic process that follows the flow between stochastic jumps. In QFT, 
it is a typical situation that H = Hq + Hi where Hq is a differential operator associated 
with continuous motion while Hi is an integral operator (often linking different particle 
numbers) associated with jumps. 

4 Global Existence of Bell's Jump Process 

In this section we deal with Bell's jump process introduced as model B in the last 
section. As we have shown in |GT03j . this process exists globally in time. In fact, for 



I Im vD t *(g \ x) (#g)- x / 2 £^ <p(x - y) %(q) 



+ 



<?t{q\x,q) 
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our proof it is not relevant whether the configuration corresponds to the fermion number 
operators. We only need that Q is any countable set, and P( - ) a PVM on Q acting on 
Jtf. In fact we can allow that P( ■ ) is a positive-operator- valued measure (POVM), a 
concept somewhat weaker than a PVM. 1 Here is our result. 

Theorem 1 Let J%? be a Hilbert space, H a self-adjoint operator on Jf? , Q a countable 
set, and P(-) a POVM on Q acting on J4f. For every initial state vector with 
l^o || = 1 satisfying 



there exists a (right- continuous) Markovian pure jump process (Qt)tm on Q with tran- 
sition rates (JTTJ) such that, for every t, Q t has distribution iz t = (&t\P( • )^t) ■ This 
process is unique in distribution. 

Some comments are necessary. First of all, how could the process fail to exist globally 
in time? Two kinds of catastrophes could occur. On the one hand, the jump rate (jllj) 
is singular at the nodes of \l/ (i.e., at such q and t for which (^> t \P(q)^f t ) = 0). While 
Qt is sitting on a configuration q it might become a node, and then the process would 
not know how to proceed. It turns out that this problem does not arise because, with 
probability one, there is no t at which Q t is a node. This is because the increase of 
the rates close to the nodes has the positive effect of forcing the process to jump away 
before the singularity time is reached. 

The second kind of possible catastrophe would be an explosion, i.e., an accumulation 
of infinitely many jumps in finite time. The main task is to show that this does not occur, 
with probability one. The standard criteria for non-explosion of pure jump processes are 
confined to transition rates that are homogeneous in time, relying heavily on the fact 
that the holding times are then exponentially distributed and independent conditionally 
on the positions; see, e.g., Sect. 2.7 of |JNor97j . This conditional independence, however, 
fails to hold in the case of time-dependent jump rates, and the singularities of Bell's 
transition rates do not allow any a priori bounds as they were used, e.g., in |Pre76| lRL53 
to exclude explosion. The only thing one knows is that the process is designed to have 
the prescribed quantum distribution (jlj) at fixed (deterministic) times, and it is this fact 
we will exploit. We will sketch our main arguments below. 

Steps towards a proof of global existence of Bell's process have also been made by 
G. Bacciagaluppi [Bac96, BD99j; his approach is, however, very different from ours. 

Concerning the technical assumptions (fT2"|) on H, P, and we note the following. 
For fixed H and P, the conditions (|T2"j) define a set of "good" initial state vectors \l/o 

1 That is, P takes values in the positive (bounded, self-adjoint) operators on 34? (instead of the 
projection operators as a PVM) and shares the countable additivity and normalization of a PVM. 



ty t £ domain(i/) 
P(q)^t e domain(if) 



\/teR, q eQ 



(12a) 
(12b) 




(12c) 
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for which we can prove global existence; this set is obviously invariant under the time 
evolution ([1)1. In fact, when if is a Hilbert-Schmidt operator (i.e., trii 2 < oo), the 
conditions (|12|) are satisfied for all POVMs P and all ^o G J^; this is also true when H 
is bounded and Q is finite. (Usually, the Hamiltonian of a lattice QFT is not Hilbert- 
Schmidt but can, at least, be approximated by Hilbert-Schmidt operators. And it is 
not unusual in quantum field theory that Hamiltonians need to be "cut off" in one way 
or another to make them treatable, or well-defined at all.) Condition (|12b|) ensures that 
the expression P(q)HP(q')ty t can be formed, and thus that ([TTjl is well-defined whenever 
q' is not a node. 

The main construction is obvious. Starting from any fixed initial time to, the process 
Qt can be expressed for t > to in terms of T k and X k , the time and the destination of 
the k-ih jump after time to, via 

Qt = Xk if Tfc < t < Tk+i 

with T := to an d X = Qt - The variables and Xk are defined by their conditional 
distributions: 



G dt, Xk+i — q\To, Xq, . . . , T k , Xk) — 

l{T fc <t} <?t (<L X k ) exp (- J a s (Q, X k ) ds^j dt, 

where the role of the "failure rate function" is played by 

cr t (Q,q) = a t(q',q) , 

g'eQ 

the total jump rate, to whatever destination, at q. 

Here is the reason why the process cannot run into a node. By definition, the 
conditional probability of remaining at q until at least time £2, given that Qt t = q, is 



rt-2 

cxi) | - / cr<( 
ti 



Q, q) dt) . 



We want to show that this probability vanishes whenever q is a node at ti but not at 
any t with t\ < t < ti- Ignoring some technical subtleties, this can be derived by the 
following simple calculation. Since a sum of positive parts exceeds the positive part of 
the sum, we conclude from (JTTJ) that 

{Q , _ V Him (y t \P(q)HP(q>)y t )]+ [lhn(V t \HPW)* t )] + 
^ ,q) £g (%\P(q')%) ~ (%\P(q')%) 

Omitting the positive part and using (JSj) we find 
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at every t with (^/ t \P(q')^/ t ) > 0. Hence, by the fundamental theorem of calculus, 



/^(ag) dt>-log<* ta |P(g)* ia )+log(* tl |P(g)* tl ) = oo, 
Jti 

since q is a node at t 2 (so that the first term is infinite) but not at t\ (so that the second 
term is finite). 

Another part of the proof we would like to sketch here is the core of the argument 
why the jump times cannot accumulate. As a convenient notation, we introduce an 
additional "cemetery" configuration oo and set Q t '■— oo for all t after the explosion 
time sup fc T fc . Let S(ti,t 2 ) be the number of jumps that the process performs in the 
time interval [ti,^]- The random variable S(tx,t2) is either a nonnegative integer or 
infinite. Our assumption ()12c|) implies in fact that S(ti,t 2 ) has finite expectation, for 
all ti < t 2 , and thus is finite almost surely. To see this we use the equation 

ES(t 1 ,t 2 )= P J2 °M)PtU)to, (13) 

with pt(q') = P(Qt = q')- Intuitively, equation (|T3^I can be understood as follows. 
a t{q, q') Pt(q') dt is the probability of a jump from q' to q in the infinitesimal time interval 
[t,t + dt\. Summing over q and q', we obtain (the expected total jump rate and thus) 
the total probability of a jump during [t,t + dt}. Integrating over t, we obtain the 
expected number of jumps. The point of equation (fT^j) is that its right-hand side involves 
exclusively the one-time quantities a t and p t . Now, one can deduce from the definition 
of the process that 

Pt{q) < (%\P(q)%), (14) 

where any case of strict inequality would have to go along with a positive probability 
f(Qt = oo) of accumulation before t. Combining (jT5j) and (JHJ), we obtain 

ft 2 



ES(t u t 2 )< f 2 J2 <Tt(q,q'){*t\P(q')%)dt 
h q,q'eQ 

E [llra(%\P(q)HP(q')%)] + dt 

q,q'£Q 

I r J] \(%\P(q)HP(q')%)\dt<oc 



q,q'£Q 

by assumption ()12cj) of Theorem 1. Indeed, this reasoning more or less dictates the 
assumption ()12c|) . 



5 Other Global Existence Questions 

Variants of the reasoning in the previous paragraph could be used in the future in other 
global existence proofs. Here is an example concerning the global existence of Bohmian 
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mechanics (with fixed number of particles). This was first proved in [BDGPZ95 under 
suitable assumptions on the potential V and the initial wavefunction One way in 
which a solution Q t of (fTUJl may fail to exist globally is by escaping to infinity (i.e., 
leaving every bounded set in M. 3N ) in finite time. To control this behavior, we suggest to 
consider an analogue of S(ti, t 2 ): Let D(t\, t 2 ) be the Euclidean distance in M? N traveled 
by Qt between t\ and t 2 . Then 

E Dfa, t 2 ) = I' dt [ dq\v t (q)\p t (q), (15) 



where v t is the Bohmian velocity vector field on M. 3N , with i-th component given by the 
right-hand side of (jl(J|) . and the expectation is taken over the randomness coming from 
the initial configuration. If this expectation can be shown to be finite, D(ti, t 2 ) must be 
almost surely finite. Using pt(q) < |\I/t(g)| 2 , the analogue of (JHJ), we see that the escape 
to infinity is almost surely excluded provided that 

dt I dq \% (q) W t (q)\ < 00 whenever t x < t 2 . 



This is a condition analogous to (|12cj) ; it is almost equivalent to the assumption A4 
of [BDGPZ95 . The proof there, however, is different, using skillful estimates of the 
probability flux across suitable surfaces in configuration space-time M 3Ar x R surrounding 
the "bad" points (nodes, infinity, points where \1/ is not different iable). The above 
argument based on (|15|) might contribute to a simpler global existence proof. 

The global existence question is still open for the Bohm-Dirac law of motion [Boh53 , 



BH93, , a version of Bohmian mechanics suitable for wavefunctions i/j obeying the Dirac 
equation. The Dirac equation is a relativistic version of the Schrodinger equation and 
reads 

i^"^" = ~ ^ {ichaj ■ Vii> + f3imiC 2 ^) (16) 
i=i 

where ipt is a function M? N — > (C 4 )®^, c is the speed of light, m.j the mass of the i- 
th particle, ct{ the vector of Dirac alpha matrices acting on the i-th. spin index of the 
wavefunction, and (3i the Dirac beta matrix acting on the i-th index. The Bohm-Dirac 
equation of motion reads 



dQ t ,i _ TptOijH 

dt ~ v>M 



(Qt,i,...,Qt,N) (17) 



where <f)*ip is the inner product in Dirac spin-space. Since these velocities are bounded 
by the speed of light, the question of escape to infinity does not arise here. But the 
question of running into a node does because, like the minimal jump rate (jHJ), the velocity 
formula ()17|) is ill-defined at nodes. 

This question can be treated in a way analogous to the previous arguments based 
on (|T3*j) and (fTo^) . To this end, let L(ti,t 2 ) be the total variation, between t\ and t 2 , of 
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log \ipt{Qt) | 2 ; i n other words, 



L{h,t 2 



t-2 



dt 



dt 



It takes values in [0, oo] and is infinite if the trajectory Q t runs into a node between ti 
and t 2 . This must be a null event if L(£i,i 2 ) has finite expectation; for the latter we 
have the formula 



EL(ti,t 2 ) 



dt 



dq 



R3iV 



0_ 

dt 



+ v t (q)- V log|^ t (g)| 2 p t (g), 



where Vt is the vector field on R 37V whose z'-th component is the right-hand side of (|17|). 
Using < (^(g)! 2 we obtain 

EL(ti,t 2 ) < f 2 dt [ dq (l + Vt ( q ).v)\Mq)f 
Jti Jr 3N Kc,t J 

As the velocities are bounded, the last expression is at most 



H 



dt 



dq 



d_ 

dt 



+ c 



Inserting (|T6|) and using that the alpha and beta matrices have norm 1 we see that this 
in turn is not larger than 



dt 



ti 



dg 2c 



^m jC 2 )^ t + 2c 



Since 



1, the last integral coincides with 



!(yVc 2 )(t 2 -ti) +4c [ dt f dq 



The question remains under which conditions on ipo the last term is finite for arbitrary 
t\ < t 2 . This is presumably the case when ipo lies in Schwartz space (containing all 
smooth functions / such that / and all its derivatives decay, at infinity, faster than 
for any n > 0). To work out the proof remains for future research |TT04| . 

While global existence of the Bohm-Dirac trajectories can presumably be proved 
with the same methods as used in BDGPZ95 for (estimating the flux across 
suitable surfaces surrounding the bad points), such a proof requires a lot of effort. It 
seems that the argument just sketched would be much easier and more elegant. 

Another global existence question that is still open is that concerning the process 
defined in |t)GTZ03l] and described above in Part C of Sect. 3, and for similar processes, 
on a configuration space that is a disjoint union of manifolds, following deterministic 
trajectories interrupted by stochastic jumps. 
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6 Deterministic Jumps and Boundaries in Configu- 
ration Space 

In this last section we describe another application of minimal jump rates, one that 
has not yet been discussed in the literature and that raises some questions for fur- 
ther research. Suppose that the configuration space Q is a Riemannian manifold with 
boundaries, or more generally the disjoint union of (at most countably many) Rieman- 
nian manifolds with boundaries. We write Q = dQU Q° where dQ denotes the boundary 
and Q° the interior. 

We develop below an analogue of Bohmian mechanics on Q, consisting of smooth 
motion interrupted by jumps from the boundary to the interior or vice versa. The jumps 
from dQ to Q° are deterministic and occur whenever the process hits the boundary. The 
jumps from Q° to dQ are stochastic, and their rates are fully determined by requiring 
that (i) the process is Markovian and equivariant, and (ii) the construction is invariant 
under time reversal, in that the processes associated with \P 4 resp. ty*_ t are reverse to 
each other, in distribution. These rates are, in fact, another instance of minimal jump 
rates. 

Configuration spaces with boundaries arise from QFT if a particular kind of "ultra- 
violet cutoff" is applied, which can be regarded as corresponding to smearing out the 
charge of an electron over a sphere rather than a ball. Here is an example. Consider 
again electrons and photons, with the electrons fixed at locations given by the finite set 
i] C M 3 , and suppose that photons cannot get closer than a fixed (small) distance 5 > 
to an electron, as they get absorbed when they reach that distance. Thus, the available 
configuration space is 

Q = {qeMW i ):d(q,T ] )>5}, (18) 

where d(q, 77) is the Euclidean distance of the finite sets q and r\. The space Q is a 
countable disjoint union of Riemannian manifolds with boundary, 



Q={J{qeQ:#q = n}. 



n=0 

Its interior is Q° = {q G T^(R 3 ) : d(q,i]) > 5}, and its boundary dQ = {q G 
For this or any other configuration space with boundaries, the law of motion 

^i = MQt) = l lm ^li m (19) 

dt y ^ ' m y*^ t y ^ ' y ' 

on Q° must be completed by specifying what should happen to the process at the time 
t when it reaches the boundary. (No specification is needed, however, for what should 
happen when two photons collide, as this has probability zero ever to occur.) The 
specification we consider here is a deterministic jump law 



Qr + = f(Qr-) 
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for a fixed mapping / : dQ — > Q° . In our example (|T5j) . the obvious choice of / is 

/(g) = { x eq: d(x,r}) > 5}, 

which means that all photons having reached the critical distance 5 to some electron 
disappear. 

Since we want the theory to be reversible, we must also allow for spontaneous jumps 
from interior points to boundary points. Since we want the process to be an equivariant 
Markov process, the rate for a jump from q' G Q° to a surface element dq C dQ must 
be, as one can derive, 

/i , [n(q).v t (q)\%(q)\*] + 
<^t{dq,q) = - \m t (q>)\ 2 n d ?>?)> ( 20 ) 

where n(q) is the inward unit normal vector to the boundary at q G dQ, the dot • denotes 
the Riemannian inner product, and z^(dg, q') is the measure-valued function defined in 
terms of the Riemannian volume measure /i on Q and the Riemannian surface area 
measure A on d Q by 

A(flnr'W)) 
" (B ' ?) Md?) ' 

with B C <9<2, and the right-hand side denoting a Radon-Nikodym derivative (the 
existence of which we presuppose). The measure z/( ■ , q') is concentrated on the subset 
f^ 1 {q') of the boundary for almost every q' . For a probability distribution on Q having 
a density function p with respect to p, one obtains the following probability transport 
equation at q' G Q°: 

?£L(qi) = _V • (p t v t ) (q') - a t (dQ, q>) p t (q') + [ u(dq, q') \n(q) ■ v t (q) p t (q)} " (21) 
cr Jag l j 

For equivariance we need that ()21|) . with |\l/i| 2 in place of p t , has the structure of the 
transport equation for \^t\ 2 that follows from the Schrodinger equation 



AI2 



d\Ml ')\ 
dt 



lm** t (q')(H* t )(q'). (22) 



This is not automatically the case, but it follows from (and therefore suggests) the 
following boundary condition relating ^t\dQ to ^|q°: for all q G dQ, 

n(q)-V%(q)= 1 (q)%(f(q)) (23) 

where 7(g) is any complex number. 2 This condition prescribes the normal derivative 
of the wavefunction on the boundary. Some boundary condition would be needed any- 
way to define the evolution of the wavefunction, i.e., to select a self-adjoint extension 



2 More generally, if ^ takes values not in C but in a higherdimensional complex vector space, 7(5) 
would be a C-linear mapping from the value space at f(q) to the value space at q. 
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of the Laplacian; whether (J2*3")l actually suffices for this, we have to leave to future re- 
search. Note that (J23j) is a linear condition and thus defines a subspace of the Hilbert 
space L 2 (Q). From and (fT9"j). one obtains equivariance with respect to the formal 
Hamiltonian H = — ir-A + Hj, where 



Furthermore, the jump rate (j2~Uj) is indeed the minimal jump rate (JBJ) associated with 
R u thanks to (1231) . 

Let us emphasize the following aspects. First, starting from the picture of a piece- 
wise deterministic process that jumps whenever it hits the boundary, we arrived with 
remarkable ease at the probability transport equation (|21jl and thus at the boundary 
condition (|23p. Second, we have derived what the interaction Hamiltonian Hi is; once 
the destination mapping / and the corresponding coefficients 7 had been selected, there 
was no further freedom. Third, it turned out that the minimal jump rate is the only 
possible jump rate in this setting. Its very minimality plays a crucial role: a jump 
to a boundary point q at which the velocity field is pointing towards the boundary, 
n(q) ■ v t (q) < 0, would not allow any continuation of the process since there is no trajec- 
tory starting from q. The problem is absent if the velocity at q is pointing away from 
the boundary, n(q) ■ v t (q) > 0. (We are leaving out the case n(q) ■ v t (q) = 0.) On the 
other hand, jumps from q to f(q) cannot occur when v t (q) is pointing away from the 
boundary since in that case there is no trajectory arriving at q. Thus, the jumps must 
be such that at each time t, one of the transitions q — > f(q) or f(q)^>q is forbidden, 
and the decision is made by the sign of n(q) ■ v t (q). 



/ /i(dg') / v(dq,q')<S>*(q')r(q)*(q) 

JQ JdQ 
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